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ABSTRACT 

This  report  presents  an  analysis  of  supercritical 
flow  in  an  overfall  whose  nappes  are  subjected  to 
different  constant  pressures.   It  contains  computations 
for  brink  depth  and  approximate  equations  for  the  nappes 
In  the  neighborhood  of  the  brink  and  far  downstream. 


ill 


1.  Introduction 

The  shape  of  the  nappes  and  the  brink  depth  are  character- 
istics of  an  overfall  which  have  intrigued  several  analysts  who 
have  used  numerical  methods  and  various  other  approximation  pro- 
cedures to  explain  the  falling  flow.   The  purpose  of  this  paper  is 
to  investigate  a  two-dimensional  overfall  under  the  assumptions, 
not  adopted  hitherto,  that  the  pressures  on  the  nappes  are  differ- 
ent but  constant  and  that  the  upstream  flow  is  supercritical. 

If  the  acceleration  due  to  gravity  is  g,  if  the  upstream 
depth  is  a  and  if  the  velocity  there  is  u,  then  the  flow  is  said 
to  be  critical  if  e  =  ga/u   is  equal  to  unity.   Our  analysis  is 
based  on  the  supposition  that  e  is  small  so  that  the  upstream 
velocity  (u  >•  >•  ^/ga)  of  the  flow  is  supercritical.   Such  a  flow 
can  prevail  in  a  region  of  space  where  g  is  small  compared  with  a 
and  u.   It  can  obviously  also  prevail  if  u  is  sufficiently  large, 
or  if  a  is  so  small  that  the  overfall  degenerates  into  a  falling 
sheet. 

We  show  below  that  the  velocity  potential  of  the  flow  must 
satisfy  a  provocative  nonlinear  integral  equation  involving  the 
parameter  e.   Taking  e  to  be  small  we  examine  the  consequences  of 
the  assumption  that  the  solution  of  the  integral  equation  possesses 
an  asymptotic  expansion  in  powers  of  e.   We  find  that  the  second 
order  approximation  provides  results  which  are  sufficiently  accu- 
rate if  e  is  small.   Even  for  the  critical  case  e  =  1  it  turns  out 
that  our  analysis  gives  a  depth  b  =  (.70)a  in  the  neighborhood  of 
the  brink  when  the,. pressure  difference  "p  for  the  nappes  is  zero. 


_^._  ^. „^_       ...  .verfall  Is  of  some  Importance  for  metering 

channel  flows  and  (.71)a  Is  the  average  oi'  approximate  values 
found  fr       application  of  other  methods  to  the  case  p  =  0. 
Besides  a  focu.  -.  :^,^...  ..n  on   brink  depth,  we  present  second 
order  equations  for  the  shapes  or  the  nappes  in  the  nelghDornood 
of  the  brink  and  far  downstream.   We  conclude  with  brief  remarks 
a.^^i  the  overfall  in  the  form  of  a  thin  sheet. 


2.  Formulation 

A  two-dimensional  waterfall  or  free  overfall  is  defined  by 
its  flow  in  a  plane  perpendicular  to  its  linear  edge.   Figure  2.1 
shows  such  a  flow  referred  to  a  rectangular  coordinate  system  whose 
origin  is  at  the  brink.   The  negative  part  of  the  x-axls  coincides 
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FIG.  2.1 


with  the  sill  or  bottom  L-,  and  the  y-axis  is  positive  upwards. 
The  domain  D  of  the  cross  sectional  flow  is  bounded  by  L-,  ;  and 
the  curves  F-,  and  To  which  define  the  free  surfaces  or  nappes  of 
the  overfall. 

The  difficult  basic  problem  is  to  deduce  the  shapes  of  r,  and 
Tp  as  consequences  of  the  assumptions  that  gravitational  attraction 
is  the  only  body  force  acting  and  that  the  flow  is  a  steady  irrota- 
tional  one  of  an  inviscid,  incompressible  fluid  whose  density  5  is 
constant.   If  v  is  the  velocity  vector,  conservation  of  mass 
requires  V-v  =  0  and  this  with  the  vorticity  condition,  Vxv_  =  0, 
insures  the  existence  of  a  function  of  the  complex  variable 
z  =  X  +iy,  namely  the  complex  velocity  potential, 

(2.1)  w  =  <|)(x,y)  +  i^{x,y)  =  F(z) 

which  is  analytic  for  z  in  D.   This  function  is  such  that  the  gra- 
dient of  its  real  part  gives  the  velocity  vector  v  =  V  Re  F(z)  =  V(j); 
and  such  that  its  imaginary  part,  ?//  =  Im  F(z),  is  the  stream  func- 
tion for  the  flow.   In  other  words,  the  horizontal  and  vertical 
velocity  components  are  given  respectively  by  the  real  and  imaginary 
parts  of 


F  '  (z)  =  i  -111/     =  Tp    +  ih 
^  ■'    ^x    ^x    ^y    ^y 


If  n  is  the  unit  outward  normal  to  the  boundary  of  D  and  if 
s  denotes  arc  length  along  the  boundary  then  the  kinematic  boundary 
condition  requires  that  the  normal  velocity 


i    =  if 


must  ce  zero  aioiig  L,,  r,  and  To*   It  follows  from  this,  by  Inte- 
gration with  respect  to  s,  that  L,  plus  r,  constitutes  a  stream 
line  which  we  can  assume  Is  defined  by  the  level  curve 

^(x,y)  =  0  . 

With  this,  the  upper  nappe  curve  Tp  is  defined  by  the  level  curve 

Tp{x,y)   =  h  =  au 

where  a  is  the  upstream  depth  of  the  flow  at  x  =  -oo ,  and  u  is  the 
velocity  there.   In  addition  to  the  kinematic  boundary  condition, 
the  energy  integral  of  the  momentum  equation 

^—  5 

5  -TTT  =  6v_- Vv_  =  -p  W'V  =  -Vp  -  g5Vy  , 

namely. 


(2.2)  ^2^^2  ^  ^  +  2gy  =  c 


must  be  satisfied  along  r,  and  r^-   In  this  equation  g  denotes 
the  gravitational  acceleration  and  p  denotes  the  pressure  in  the 
fluid.   Since  the  flow  is  supposed  to  be  steady  the  constant  c 
remains  the  same  no  matter  where  the  left  hand  side  is  evaluated. 
Therefore  with 


X 


(-00  ,y)  =  u 


't'yC-oo  ,y)  =  0 


p{-oo,a)  =  pg 


we  have 


p        2pp 
(2.5)  c  =  u''  +2ga  +  -^ 


Under  certain  circumstances  the  pressure  Pp  acting  on  the 
upper  nappe  of  the  overfall  is  different  from  the  pressure  p, 
acting  on  the  lower  nappe.   For  example,  wind  action  on  the  upper 
nappe  or  vapor  pressure  action  under  the  fall  could  cause  pressure 
differences.   The  pressures  Pp  and  p-,  are  not  likely  to  be  the 
same  if  the  overfall  is  confined  between  parallel  walls  perpendicu- 
lar to  the  brink  and  the  fall  is  over  a  vertical  cliff  so  that  the 
air  under  the  fall  is  isolated  from  the  rest  of  the  bounding  air. 
For  a  case  like  this,  intentional  pressure  differences  could  be 
produced.   In  order  to  approximately  cover  these  possibilities  we 
suppose  that  the  prescribed  pressures  Po  and  p-,  are  constants  such 
that  Pp  -  P-,  is  not  large. 

Although  the  domain  D  is  not  known  at  the  outset  we  can  work 
with  a  fixed  domain  in  the  w  =  ^  +  lip   plane  if  we  introduce  the 
inverse  function 

(2.4)  z  =  x{^,f)   +iy((t),^)  =  G(w)  . 

Under  the  mapping  from  the  z -plane  to  the  w-plane,  D  is  mapped  into 
the  strip  S  defined  by 


0  <  ■'    "-  w  <   h 


XI  <  ^   =  Re   vi  <  oo    , 


the  or!     r  the  w-plane  to  be  the  Image  of 
the  origin  of  the  z-pla..c.  .  jln  fixed  In  this  way  the 

Images  of  L,,  P-,  and  r^  are  respectively  '"  .  ''-,  and  Zo  as  shown 
in  Fig.  2.2.      :e  there  are  no  stagn      points  in  the  overfall. 
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the  function  G{w)  must  be  analytic  in  S  and  we  have 


FIG.   2.2 


dz 
3w 


G'(w)    =   x^  +  iy^  =    -i^^  +  y^ 


—  =  1  ^  1 

dw       "5    +  lip         ^    -  i<t) 

^     ^      ^      ^      y 


^x        ^y 


from  whi ch 


G'{w)G'{w)    =    |G'(w)| 


2   .     2 


=  '^^  '='*  =  FTP 


and 


^x     ^y 


Re  G ' ( X ) 


G'(w) 


;    i 


Im  G' (w) 


^x    i.,,..m"?  '    ^y  ~  i..,_.m2  • 


G'(w) 


These  relations  with  (2.2)  and  (2.3)  show  that  G(w)  must  satisfy 
the  following  boundary  conditions: 


G'((t)) 


Im  G'  ((f))    =   0    , 


1                                    ,              P                    2(pp-p.  ) 
^^ — o   +  2g   Im  G((t))    =  u""  +2ga   +  ^ — — 


^  +   2g   Im  G((t)  +  ih)    =  u^  +  2ga    , 


G'  (i)  +ih)  I 


(f)    <   0    ; 


<})   >■   0    ; 


-CD      <    ({)     <     00 


We  can  now  map  into  any  domain  we  think  may  be  helpful  for 
the  analysis,  whether  strictly  numerical  or  not.   Instead  of  the 
strip  S  in  the  w-plane  we  are  going  to  transform  the  problem  into 
one  for  the  upper  half  of  the  complex  C-plane  into  which  S  is  mapped 
by 


C  =  -e 


-TTw/h 


Under  this  mapping  the  images  of  M-,  ,  Z-,  ,  and  2p  are  respectively 
L^,  r,  and  Fp  as  shown  in  Fig.  2.3.   The  image  of  the  brink  0  in 


L' 


-I 


FIG.  2.3 
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the  z-plap.e  Is     .t  C  =  -1  while  the  source  and  sink  of  the 

c.     X  .c'j^:;e  respectively  the  point  at  Infinity  and  the  origin 


in  the  C-plane.  We  also  have 


G(w)  =  G[^  in  (-  ^)]  =  0(C) 


G'(w)  = 


-  ^o'(a  . 


The  problem  now  Is   to   find   the   function  n(C)   which   Is  analytic   for 
"      '   >   0  and   satisfies   the  boundary  conditions 


(2.5) 


ImO'(e)    =   0    , 


^  <  -1  ; 


(2.6) 


h  P  — 

g   g,     . ?  =   u''+  2ga(l+  p)  -  2g   Im  0(U    ,         -1   <    ^    <   0    ; 

7r^rio'(e)r 


(2.7) 


^   A  ^ Tj  =  u^  +  2ga  -  2g  Im  0(  ^ )    , 


0   <   ^    ; 


where 


P  = 


P2-P1 
5ga      • 


For   the   neighborhood  of   Infinity  -ttCO'CO    must  behave   like 


-TrCn'(C)    ~  hG'(w) 


^=- 


h(^^       +       i<t>y) 

2  .   i2 


CO 


^x     ^y 


x=-oo 


^  ^  =  a 


By  integrating  the  condition  ImO'(^)  =  0  we  find  ImO(^)  =  const, 


8 


for  -00  ^    i   '^    -1,    but  since  this  part  of  the  real  axis  corresponds 
to  the  sill  where  y  =  Im  G{^)    =   Im  0(4)  is  zero  the  constant  of 
integration  is  zero  so  that 

Im  0(U  =  0  ,  -co  <  ^  <  -1  . 

This  implies  that  0(C)  can  be  analytically  continued  across 

Im  C  =  0  ,  -oo  <  Re  C  ^  -1  ; 

in  accordance  with  the  formula 


o{C)  =  o(C)  . 

Therefore  0(C)  must  be  analytic  in  the  cut  plane 

A:   0  <  arg  (C+1)  <   2tt 

subject  to  a  nonlinear  boundary  condition  along  the  banks  of  the 
cut. 

In  the  beginning  the  problem  was  to  determine  a  function  F(z) 
analytic  in  an  unknown  domain  D  but  subject  to  two  boundary  condi- 
tions along  each  of  the  unknown  parts,  T-,  and  Tp,  of  the  boundary 
of  D.   We  have  seen  how  the  dual  boundary  conditions  enable  us  to 
transform  the  problem  into  one  for  the  determination  of  a  function 
0(0  analytic  in  a  known  domain  but  subject  to  one  nonlinear 
boundary  condition.   The  nature  of  the  nonlinearity  is  such  that 
the  probability  of  finding  n(C)  in  closed  form  is  practically  nil. 


Aside  from  a  strictly  numerical  analysis  It  seems  that  the  best 
t;     :n  be  done  Is  equivalent  to  expressing  n(C)  as  a  function 

h  satisfies  a  nonlinear  Integral  equation  and  then  resorting 
to  some  approximation  scheme.   Various  Integral  equations  can  be 
derived  depending  on  the  part  of  the  nonlinear  boundary  condition 
we  think  of  as  knovm.   We  proceed  to  regard  the  right  hand  sides 
of  (2.6)  and  (2.7)  as  knovm  and  derive  an  Integral  equation  for 

in[-7rCO'(C)]. 

Consider  the  function  £n["'^^^'  ^^^]//Z+l.      This  function  is 
analytic  in  A;    and  along  -oo  <  ^  <  -1  we  have 


J_  arg[^lSV(lil 


J"U^        " 


=  0 


because  Im  O'  (^)  =  0  for  -oo  <  ^  <  -1  and  -win'  (^)  —  a  as  i  -*■  -co 
Hence  by  continuation 


1    ^^r-7r(^-lo)n'(e-io).  ^      1    ^^^-7r(e+lo)0'(^+lo)^ 


Now  an  application  of  the  Cauchy  integral  formula  to  the  slit  plane 
A  gives 


10 


^^^-7r(A+lo)Q'(A+lo)-|       ^ 


00 

1  _^_^^ 

h        ■■        27ri 


■    + 


j,^[-Tr(A+io)n'(A+lo)-j 


or 


00    „      „    r-7r{A+lo)0' (  A+1o)t 

*"  ""^        '{^  /ATI    (A-C) 

Then  from  the  conditions  (2.6)  and  (2.7)  we  have 


°  2 


i/C+1       r       ^n[u^  +  2ga(l+?)  -  2g   Im  n(A)]    ^^ 
^"^       ^1  yATl    (A-a 

(2.8)       ^n[-^^g'^^h    = 

^ 00       „ 

^  i/Z+1      r       in[u^  +  2ga  -  2g  Im  n(A)]  ^^ 


0 


/A+l  (A-C) 


In  the  above  we  have  implicitly  assumed  that  in[ -ttCO' ( C)  ]  possesses 
no  singularity  at  C  =  -1.   If  we  introduce 


^  "   2 
u 


and  use 


oo 


l/C+1   r       dA     ^  _  1 
^^   -(_^  /A+T  (A-C)    ~  ^ 


then  (2.8)  becomes 


11 


0    .    r-,    .   o./i  .-\       2e 


in[l  +  2e(l+p)  --1^  Im  0(A)] 
"""W    J 

(2.9,     in[^imC)u^  ^ 


h 

,/rTT      r°°    in[l+2e-ii  ImO(X)] 

dA 


^^/ 


/mT  (a-C) 


Finally,  if  we  substitute 

h  =  au  , 
0(0  =  ai  +  aa)(0  , 

and  use  the  unit  function 

/I     t  >  0 
H(t)  = 

[  0      t  <  0  ; 

we  see  that  (2.9)  can  be  written  as 

(2.10)   £n[-7rCa.'(0]  =i^  f       ^n(l  +  2e[p^  ( -A)  -  Im  co(  A)  ]  ]  ^^  ^ 

An  integral  equation  for  a>(0,  -1  "^  4  "=  oo  ,  can  now  be  produced  by- 
letting  C,  -*   ^+io.   However,  we  do  not  need  this  equation  explicitly 
because  we  are  interested  in  studying  the  consequences  of  the 
simple  and  seemingly  most  natural  assumption  that  a)(C)  possesses 
the  asymptotic  expansion 

(X> 

03(0  =  ZZ  eV(C) 
for  small  e  and  at  least  \l^\    >■  r  >  0. 
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3.  Terms  of  the  Expansion 


The  substitution  of 


(3.1) 


00 


0(C)  =  ai  +  acD(C)  =  ai  +  a  ^^  e'^a)^(0 


n=0 


in  the  boundary  conditions  (2.5),  (2.6)  and  (2.7)  would  generate 
an  infinite  sequence  of  boundary  value  problems  which  could  be 
solved  recursively.   This  sequence  is  equivalent  to  what  is  found 
by  substituting 


(3.2) 


CO 


n=0    ^ 


in  (2.10)  and  equating  coefficients  of  like  powers  of  e.   The 
result  of  doing  this  gives  the  following  for  cd  ,  cd-,  and  oOp: 


(3.3) 


^n[-TrCcD;^(a]  =  0 


(3.4) 


"^1^^^    i/^+T   r°°  [PM-(-^)  -Im  cOq(A)] 


%t?) 


77" 


dA 


-1 


/A+1  (A-a 


(3.5) 


co'(a 

o 


{Im  CD^(A)+  [pM.(-A)  -  Im  a3^(A)]^} 
/ATI  (A-C) 


dA 


It  should  be  noted  that  since  0(C)  is  analytic  in  the  slit  plane  A, 
this  must  also  be  the  case  for  each  en  (C).   Moreover,  because  the 
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luilvatlve  a)'(s)  appears  iw  li.e  Integral  equation  (2.10)  we  need 

"     (O-   Such  Is  provided  by  our  mapping 
i'hloh  requires  C  =  -1  to  correspond  to  2  =  0.   Hence 


00 
0  =0(-l)  =  ai+a  J~   e"a)  (-1) 


fr  —  '.-hlch 


CO 


co^(-l)  =  -i 


(D^(-l)  =  0  ,     n  7^  0  . 


Integration  of  (3' 3)  yields 


(3.7)  0)^(0  =  -  ^  £n  C 


so  that  along  the  real  axis 

-1     -1  <  ^  <  0 
0      0  <  ^   . 


(3.8)  Im  (o^C^)  = 


The  function  to  (C)  corresponds  to  the  sill  based  Jet  that  could 
exist  if  no  gravitational  force  were  acting.  If  g  were  zero  we 
would  have 


(3.9)  z  =  0(0  =  ai  +  acD^(a 


=  ai  -  -  £n  C  . 


For  -1  <  ^  ■<  0,  corresponding  to  the  lower  nappe,  (3-9)  gives 


14 


(3.10)  y  =  Im  z  =  Im  0(?)  =  a  +  a  Im  a)^(  4)  =0  ; 

and  for  the  upper  nappe,  corresponding  to  0  ^  4  "^  oo  j  (3'9)  gives 

(3.11)  y  =  Im  z  =  Im  O(^)  =  a  . 
the  velocity  components  of  the  flow  (3.9)  Q-re 


(3.12) 


^  ^  Re  G'(w)  ^  _^  Re  TrCo'(a  ^  h  .  ^ 
^    |G'(w)r       |7rCO'(C)r    ^ 


^  .  _h  ^^-^^'iO   =  0  . 
^      |TrC0'(C)r 


These  results  describe  a  uniform  horizontal  flow  with  velocity  u. 
We  are  interested  in  the  perturbation  of  this  flow  which  is  caused 
by  gravitational  action. 

Expansions  different  from  (3.1)  can  be  assumed;  for  example 
an  expansion  parameter  different  from  e  -   ga/u  is  a  possibility. 
Besides,  we  can  contemplate  Iteration  processes  for  the  determina- 
tion of  OD.   However,  in  our  expansion,  oo-,  ,  as  given  by  (3.4), 
measures  the  first  order  effect  of  the  presence  of  the  gravitational 
acceleration  g.   The  substitution  of  (3.7)  in  (3.^^^)  yields 

a,'(C)  =  -liP±llZ^  r  dA 

^  ir'^fi  ^_^   /ATT  (A-C) 

and  after  performing  the  integration 

(3.13)         a.;(C)  =  -  iiftll  ^n /'^^SiLlI 
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__;.je  we  must  have  a>,(-l)  =  0,  (5'13)  gives 
(5.14)      a^(0  =  -  ^^  /  ^  ^n(;^7Y) 


dX 


where  the  path  of  integration  from  -1  to  C  is  in  the  slit  plane  A. 
As  far  as  we  know  the  integral  in  (3.14)  is  neither  tabulated  nor 
expressible  as  a  finite  combination  of  tabulated  functions.   In 
order  to  study  ai,(C)  in  the  neighborhood  of  various  values  of  C, 
without  numerical  integration  it  appears  that  we  must  resort  to 
several  different  expansions  of  (3.14). 

The  integral  in  (3.l4)  can  be  written 

r  1  in(  /^^~  I   )dA  =  f    i  [in  A  -2  in(v/A+r  +  1)  ]dA 

C  0 

=  /  i  (£n  A-2  £n  2)dX  -  2  f     [  ^n(/ATr+ 1)  -  £n  2]  ^^ 
-1  -1 

,c         

-  2/  [in(yA+l+^l)-  £n   2]  ^^  ^ 
0 

the  value  of  the  definite  integral  is 

1 


_2  f     [-g"(/>^+  1)  -  ^n  2]  d^  ^  r 


in  A  dA 


_1  0  ^^"^^-^  {/TTa  +  D 

2 


=  2  £n^2  -  ^  . 


For  |a|  <  1  we  have 
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2    ^n(/l+A  +1)   -  2    in   2   = 


^       00     {-l)^r(|+k)         ^k 
/tt  k=l 


It   follows   that   for    |^|    *=   1   we   can  express  od,  (C)    in  the   form 


(3.15)      03^(0    =   - 


i(p+l) 
"2 — 


IT 


^   (^n    C    -    2    in   2)^   +  ^ 

CD     {-l)kr(i  +  k) 
+   27ri    in  2  +  ^  ^^  — 
A  k=l 


kJ 


^k 


If  C  is  In  the  neighborhood  of  the  origin  which  corresponds  to  the 
flow  far  downstream  from  the  brink  of  the  overfall,  od^(C)  behaves 


like 


(3.16)      co^(C)    = 


ii^i^   [|(in    C  -2    in   2)^+  ^  +  27ri    in   2]  +o(|C|  ) 


V 


3 


The  representation  (3.14)  is  equivalent  to 


(3.17) 


03^^(0 


i(p+l) 

TT 


r  CO 


-1 


lin, 


/A+1 -  1 

yA+1  +1 


dA 


00 


c 


1  in 


yA+i  -1 


dA 


The  path  of  integration  for  the  definite  integral  is  in  the  slit 
A-plane  from  -1  to  oo  —  it  does  not  pass  through  the  origin.  An 
evaluation  shows 


03 


1  ,,/yA+T- 1 

/a+I  +  1 


A 


£n 


dA  =  27ri  in  2 
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.^ ,  as  we  can  see  from  (5.1?)  that 


0^(00)  =^^i^  in   2 


For   |X|    >■  1  we  have   the  expansion 


VxA+l+iy  /ir  k=0  T      k  +  -^ 

^  ^  kJ(k  +  ^)X 


Therefore  the  expansion 

l{p+l) 


(3.18)  03^(0  =  -  ^^ 

IT 


,   OD  (-i)^r(^+k) 

27ri  in  2  +  —  7 


A  ^^^  k.'(k  +  |)^ 


C 


-i 


Is  valid  for  |C|  >■  1- 

Using  the  substitution  yX+T  =  u,  the  representation  (3-l4) 
can  be  changed  to 


/^+i 


<„  (C)  =   -  iip^ 


2u   .   #u-l^, 
-^  in  {^^)du 

u  -1 


We  find 


which  shows  that  for  the  neighborhood  of  C  =  -1>  ^n(0  can  be 

1/2 
expanded  In  powers  of  (^+1)  '  . 

l(p+l) 


(3.19)  0)^(0  =  -  ^-^ 

TT 


00  /  n 


1.  ta(-o  +  4i;^(^5ix]^^ 


2n+l 


where  the  expansion  of  £n(-0  has  been  suppressed. 
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A  power  series  expansion  for  the  neighborhood  of  ^  =  1  can 
be  deduced  from  (3. 17).   It  turns  out  to  be 


(3.20)      CD^(C)  = 


l(p+l) 


TT 


27rl  ^n  2  +  3.812 

-  1.763(C-1)  +  1.235(^-1)' 

-  .779(C-l)^  ... 


4.  First  Order  Approximation 

Before  proceeding  any  further  It  Is  of  Interest  to  examine  In 
more  detail  the  Implications  of  the  first  order  approximation 

(4.1)     z  =0(C)  ~0^(0  =  al +acD^(0 +aeoo^(z)  =  z-^tC)  • 


This  gives  the  first  order  flow  parametrically  in  terms  of  the 
complex  parameter  C: 


-i^—  =  Re  (  -  ^  in  C  +  ea)^(C)} 


yi(C) 


1  +  Im  {  -  :i  in  C  +  eco-^(0}  . 


The  first  order  shape  of  the  lower  nappe  in  the  neighborhood 
of  the  brink  is  given  by  (3. 19)  if  we  take  -1   <    11,  =    ^   <   <   0.      We 
have 
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x^(U 


yA'.) 


a 


1 


^SiSliltn   \i\ 


ii£±i).4n:fiz 


TT 


^e(p-^l) 
5~ 

IT 


n=l 


k=l 


1   \  (1+e) 


2n+l 
~"5~ 


2n+l 


'mf^   .  ^  (1.^,5/2  , 


Our  development  is  based  on  the  assumption  that  e  is  small.   Notice 
that  if  at  this  point  we  try  to  relax  the  assumption  and  suppose 
that  the  last  result  is  valid  for  e  =  1,  p  =  0,  then  we  are  left 
with  X, (^)/a  =  0.   This  is  one  more  example  of  the  complications 
which  usually  arise  when  we  investigate  critical  flows. 

The  approximate  shapes  of  the  first  order  shapes  of  the  nappes 
far  downstream  from  the  brink  come  from  (3'l6).  For  the  lower  nappe 
we  have 


x^(^) 


in    I  el     ,    £(p+l)    -en    Ul 


TT 


IT 


yi(l) 


£(P+1) 

— 5 — 

TT 


^   {£n    \i\    -    2   £n   2)^ 


TT 

"5" 


and  for  the  upper  nappe  we  have 


x-^ii) 


_    In   i   _^  e(p+l)  2  in  2 

V  IT 


yi^^^  _  ,   e(p->-l) 


TT 


I  (£n  e  -  2  in  2)^  +  ^ 
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These  are  parametric  equations  for  the  respective  parabolas 

o 

~          aE{p+l)                  e{p+l)  r~  2a    r    A—  n ^      .i    „      ^) 

y-,    =  — '-f '- ^ '- ^  •     X      -  — -  [e(p+l    -  1]    £n   2] 

^  6      2a[e(p+l)  -  1]^   *-  ■"■    ^  ^ 

and. 

As  the  same  coordinate  y-,  decreases  to  -oo  the  magnitude  of  the 
difference  of  the  abscissas  of  the  parabolas  Increases,  that  Is, 
they  diverge.   Since  this  Is  contrary  to  the  observed  behavior  of 
the  nappes  of  a  waterfall,  the  parabolas  should  not  be  accepted  as 
good  approximations  particularly  when  e  =  ga/u  is  not  small. 
However,  if  we  take  the  case  for  which  the  pressures  on  the  nappes 
are  the  same  (p  =  0)  and  let  a  ->•  0  we  find  that  each  parabola 
reduces  to 

2u 

which  gives  the  shape  of  a  thin  heavy  sheet  of  fluid  projected 
horizontally  from  the  brink  with  velocity  u  and  subject  to  equal 
pressures  on  Its  faces. 

It  can  be  seen  from  (3.1?)  that  for  the  upper  nappe,  where 
0  <  i, 


x-l(?) 


1  in  ^  +  -^lP±iI  2  £n   2 

TT  -IT 


and 
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iilii.i.ii^  /  -  Y/iii^idx 


tt"        -L  V  sA+i  + 1 


The  first  or  :  Lrlnk  depth  to  depth  upstream  at  oo  is 

give:  '  ■ 

(4.2)      '-1^   =  1  .  ^iliil  r    ^  inf^/I^L^   dX 


4        v/^^ 

^O 


where 


£n  ^   =  e(p+l)  2    i.n   2    . 
o 

We  cannot  expect  the  formula  (4.2)  to  be  very  accurate  if  e  =  ga/u 

Z't   small.   However,  it  is  interesting  to  see  what  (4.2)  gives 
If  p  =  0  and  the  velocity  far  upstream  is  critical  so  that  e  =  1. 
For  these  conditions,  which  prevail  in  many  engineering  analyses  of 
transition  from  channel  flow  to  overfall,  our  first  order  approxima- 
tion gives 


«o  =  * 


and 


n — ? —  *  T 

TT^ /?  k=0  k.'(k+i)       k  +  -^ 


4 


-I       ^^  ( o       ^^  ^   3 


;?  ^"  -  3^  "  1^00  • 


A.  .80 
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This  is  to  be  compared  with  the  value 


^1 


which  is  the  average  of  contraction  coefficient  values 

.715  Rouse  [1] 

.705  Southwell  and  Vaisey  [2] 

.720  Jaeger  [3] 

•720  Roy  [4] 

.710  Fraser  [5] 

.676  Hay  and  Markland  [6] 

deduced  from  various  experimental  and  mathematical  approximation 
procedures  different  from  the  one  we  are  using. 

Physical  consideration  of  the  decreasing  slope  of  the  upper 
nappe  streamline  of  an  overfall  from  0  to  -00  ;  and  consideration  of 
our  choice  of  mapping  to  the  complex  w-plane,  where  (|)  =  0  corre- 
sponds to  the  brink,  together  show  that  the  point  on  the  upper  nappe 
directly  above  the  brink  must  also  be  a  ooint  on  the  potential 
curve 

H^,y)   =  -7  ,  7  >  0  . 

This  point  corresponds  to  a  point  C  =  4  ^  1  because 

TTW        TT  /    ,  . ,  \ 

c,  =    -e  =  -e 

Try_ 
=  e    >  1  . 
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*-1tut,  -   ^  1  In  (U.  2)  we  find 


y^d) 


1 


=  1  -  ^  ['^ 


2.3     5 
^   5^  "  19^ 


..] 


~   .62 

a  value  which  is  presumably  too  low.   It  is  amusing  to  notice  that 
the  mean  of  this  value  .62  and  .80  is  .71. 


5.  Second  Order  Term 

We  turn  now  to  a  study  of  the  second  order  term  cDp(C)  of  the 
expansion  (3.1).   This  term  is  defined  by  (5-5),  namely. 


(5.1) 


1 

2 


'(^{{O 


-2 


.%":). 


00 


{Im  oa-L(A)  +  [pvi(-A)  -  Im  cOq(A)]  } 


-1 


yATl  (A-C) 


dX 


Substitution  of  the  now  known  ^^^(C)  and  oo^tC)  gives 


(5.2)  CD^{0 


(p+1)^  ^ri^/'  /C^-1 

A+1+1 


27r^C 


(P  +  1)CD-[(C) 


00 


V    C 


Im  oi,  ( A)dA 


'_^  yx+i  (A-o 


which  holds  for  Im  C  >■  0.   Since  Im  co^(U  =  0  for  i   <   -1;    the 


2h 


function  co-,(C)  is  analytic  In  A,  the  slit  C-plane,  and  such  that 


cci^(^-io)  =  oo-|^{^+io)  for  -1  <  ^  <  oo  .   Also 
(5.3)  c.^(oo)  -   ^(F-^lMng  , 

as  given  In  Section  3,  is  real.   Hence  equation  (5.2)  can  be  written 


"^  27r^C  V/C+l+l/  ^ 


27r^C      _i  /A+To+T   (A-C) 


00 


[cD-,  (A-io)  -  CD-,  (oo  )] 

dA 


27r''C     ^_i  yA-io+1    (A-C) 

where  the  paths  in  the  integrals  are  along  the  upper  and  lower 
banks  of  the^cut  along  the  real  axis  of  the  C-plane  from  -1  to  +oo 
A  deformation  of  these  paths  into  the  remainder  of  the  real  axis 
yields 

(5.4)   cn'(C)  -  i^±^  ^n^f/^g-M-  (?+l)cD'(a 


-00 


+  ^ +  — ^—  /     dA 


TT^C  ^\  /ATI  (A-C) 


Integration  of  (5'^)  gives 
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r'      Loi,  (o)  -03,  (oo  )  ] 
(p-l)[cB^(C)  -co^{aD  )]  -  :^  j   -^ ^ do 


1  in  /m-^  -  ^/^ 

_^  -    .-^^  [a3i(A)-co-^(oo)]  1       /C+1  +  1 


TT^    _^  yx+i 


+  ^  ., arc  tan  -^^-^ — 


dX 


_;.  ^:-:er  to  go  on  with  the  exposition  of  a)p(C)  we  need  the 
values  of  some  definite  Integrals.   These  Integrals  are  listed  in 
the  appendix  with  indications  of  methods  for  their  evaluation. 

Recalling  that  a)p(-l)  =  a)-,{-l)  =  0,  we  see  that  for  C,   =  -1, 
(5-5)  becomes 

(oo)=iP±H!      r    iin2/2i±LLi^da   -    (p+l)a.,(a)) 


0^2 


-U 


[o),  (A)  -  CD,  ( 00  )] 

— dA 

A/A+1 


A   reference   to   the  definite   integrals   tabulated   in   the  appendix 
shows   that   the   value   of  002(0   at   infinity  in  the   A-plane   is 

o  2  

(5.6)     cDp(oo  )    =   -  iP±2J_  5  £n  2  +  iP±I)-  (8.416)    +  -^^±U   (8.4l6) 
"^  ^  27r^  TT^ 
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The  shape  of  the  upper  nappe  depends  on  the  behavior  of 
CDgCC)  for  C  =  4  ^  0.   From  (5-5)  we  find  that  the  real  part  of 
u^^i^)    for  ^  >■  0  is 


(5.7)     Re  Gi^{^) 


00 


1  ,n^U^^^^\aa 
^     V/^+1  +  1 


00 


Im  03-,  ( o)da 


+  032(00  )  ■ 


^ 


The  point  1=1  corresponds  to  a  point  in  the  neighborhood  of  (0,b) 
where  b  is  the  brink  depth  in  the  physical  plane  of  flow.   For  the 
neighborhood  of  ^  =  1  the  Taylor  expansion  of  (5-7)  is 


(5.8)   Re  o^^i^)    =   032(05) 


27r^ 


3.512  -  3-108(4-1)  +  2.800(^-1) 


1.327(4-1)^  ... 


2  n 


(p+1) 


TT 


3 


7.869  -  3.812(4-1)  +  2.788(4-1) 


_  -   2.270(4-1)^ 


2  -1 


The  imaginary  part  of  032(4)  for  4  ^  0  is 
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(5.9) 


I-U)o(4) 


(p+i) 


TT 


2  £n^2  +  ^ 


OD 


Useful  upper  and  lower  bounds  can  be  found  for  Im  (o^CO  if  we  use 
the  inequali"-  '/ -7)    in  the  appendix.   If  (p+l)  >  0  we  find 

(5.10)        Im  0)2(0    <   -(P+l)    Im  a>,{0    ^    ^P^^^f'^^^^    inf  ^^^ 

-^^^^^4^n(v/lT^+/l) 


TT^  /f 


and 


(5.11)        Im  03.(0   >   -(P+l)    lmo),(0    +  IPliHMl^l  ^^/HLj. 

^  ^  TT^  \/T+r  +  i 


l6iP+ll  J,  ^^  (^1^  +  yi)  _  .  00597 (p+l) 


The  shape  of  the  lower  nappe  in  the  neighborhood  of  the  brink 
depends  on  the  behavior  of  03p(C)  for  C  =  ^  near  ^  =  -1.   An  inte- 
gration of  (5.4)  shows  that  for  -1  <  ^  <  0 
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0)2(0 


2v^ 


i 


1  ^^2^y^Ti- 1^  ^^  _  (-+i)^^(^) 


+ 


IT 


do 


-1 


fl  ta^i.z/i±i' 


-00 


+ 


TT 


-1 


[CD-J_(A)   -  03^(00  )] 


1+/I+I 


dA 


+ 


2/-A-1  ^      .     y^+1 

— i — ^ ar>p     ran   i— 2 


arc    tan 


/^ATT  > 


The   real   part   of  ocig(^)    for   -1   <    ^   <   0   is 


(5.12)      Re  c^^U)   =    -  liPJll!  in    (-1)    +  iP±li!      , 


27r' 


-1 


,  (P+1)    f  !ij~(Y^^\io±i) 


1  +/0+1 

^2n+l)/2 


da 


ir^       ^{^    ^  n=l  V  k=l 


2k- 1  2n+l 


da 


which   can  be   expanded   into 


(5.13)   Re  0)2(4)  = 


3(P+1)' 
27r 


^n  (-U 


2{p+lf 


IT' 


li^.  5(^,1)3,  49(^,,)4 


90' 


4(p+l) 


v 


^(4.1)5/2.4(4.1)^/2.^(4.1)9/2 


9^' 


valid  for  -1  <  4  <  0. 

The  imaginary  part  of  cd2(4)  for  -1  <  4  ^  0  is 
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(5.1^)   ImojgCU  =  -  2(p+l)  Im  ui^(a 


+  iP±ll  In    {-i)    [£n  i-i)    -k    In  2] 
+  il±ll  (8.416)  Jl^J^^ 


From  the  Inequality  (A.?)  we  find  that  if  (?+l)  >   0  then 


(5.15)  ImcopC^)    <  I(^)    -  iiiPlU^  in  ( /TTf  + /?:) 

+  IP±ll   [2    in^S+^l 

TT 

and 

(5.16)  Imo)  (^)   >  I(U    -   l6(p+ll_l   ^^  (y^^  ^^) 

+  iP±|l  [2  £n22  +  I^]    -    {. 00597)  (p+1)    ; 

TT 

where 

I(^)    =   -   2(p+l)    Im  m^(^) 


+  iP±ll  £n  (-e)  [in  (-^)  -4  in  2] 


? 


27r 


+  lP±ll  (8.416)  inf  i^^^i±l 
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If  we  are  interested  In  studying  the  shapes  of  the  nappes 
far  from  the  brink  we  need  to  know  the  behavior  of  cOp(C)  as  1^  -*■  0. 
By  using  (5.5)  and  the  integrals  listed  in  the  appendix  we  find 

,5.17)     ,  ,c)  =  -  i£ii^  3  to  2  +  'P^^)''8-^6)   +  <P^1"8.1H6) 

_  iP±ll-[-  ^(^n   C   -   2   in   2)^  +  4.8o8] 
2Tr^  ^ 

+  l(I±l)^[(£n    C-2    in   2)2   +   2^] 
2Tr  ^ 

i^iil[27r2(ini;  -2  in  2)  +  ( in  C  -  2  in  2)^-12(1.202057)] 


+ 


6^3 


i(p+l)(8.4l6)(^^  ^  -  2  in  2)  -  iiP±^  +  o(UI) 


For  the  lower  nappe  with  p  =  0  (5-17)  reduces  to 


,>,      3  -^n  2    .600 

TT 

+  -i-[(in    Ul     -    2    in   2)^   -    5Tr2(in    U|     -    2    in   2)] 
37r^ 

+   i[       1   +  Ml6(^n    Ui-    2    in  2)  +  ^(in    U|  -  2    in2)2]  +o(?) 
5  TT^  2^"^ 


It  will  be  recalled  that  for  p  =  0,  and  £  =  1,  x-j^(^)  vanishes  for 
the  lower  nappe.   This  is  not  the  case  here  for  Re  (i:^^^^^ ' 
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nk  Depth 


We  proceed  now  ti. 


what    • 


: aer  approximation 


(6.1)      2   =  0{^)   ~  r.o(C)   =  al  (O  +aea)^(O  +  ae^C02{O    =   z^{i:) 

Implies  with  resp-  nk  depth  b  of  the  overfall.   The 

^. .„^.:...  .^^...  .. ^,  _  -^   of  considerable  interest  because  it 

ars  in  var"        .ulae  for  flow  measurement  when  the  overfall 
1.       .;  a  meterit.       e.   When  the  overfall  exists  in  a 
,..-..nel  the  upstream  flow  ;^  ..'^-arly  critical  and  hence  the  contrac- 
ts     • ficient  for  critical  flow  is  particularly  interesting. 

Our  previous  results  for  oo  ,  co-,  and  oo^  show  that  the  para- 
...u  ,xic  equations  for  the  upper  nappe  corresponding  to  the  neighbor- 
■  ^  =  1  are 


(6.2) 


x^(^) 


=  Re  {o^^{^)+  Eu^^ii)  +e  ca^ii)] 


1  in  ^  +  -^iP±^  2  ^n  2 

TT  IT 


+   e 


(p+1)^  ^  ;n  ?  +  (P+l)^(8.4l6)  .  (p+l)(8.4l6) 
?    ;>  ijn  ^  +      — ^      +       ^ 


£^(P+1)^ 
27r^ 


3.512  -  5.108(4-1)  +  2.800(4-1)' 
-  1.527(4-1)^  ... 


e^(p+l) 


TT' 


7.869  -  3.812(4-1)  +  2.788(4-1)' 

_  -  2.270(4-1)^  . .. 
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and 


(6.5)     -\—  =   1   +  Im  {cDq(4)    +   eto^(^)    +   e%(e)) 


1   .IP±ll[e-s2(p+i)] 


TT 


5.812   -   1.763(^-1) 


+  1.255(^-1)^  -    .779(?-l)^- 


2,- 


+ 


e^{?+l)(8.4i6)   ^^//lTT-1 
7^  V  A+l +1 


2-,^,- 


l-i^iP+il  ^n  (yiTf  +/?)    -    (.  00597  )e^(p+l)e 


tt"^/? 


where   0   is   real    such   that   0   <   0   <   1   if    (p+1)    >•   0.      In   this   section 
we   confine   ourselves   to  a  discussion  of    the   case   for  which  p  =  0. 
For  this    case   the  upper  nappe   in   the   neighborhood   of   |   =   1   is  given 

by 


(6.4) 


Xo(?) 


^n   I   ^   e2  ^n   2   ^    ^2 
+  +   e 


IT 


IT 


3   ^n   2   ^   12.624 
^3 


■w 


9.625  -   5.364(^-1)   +  4.188(^-1) 
^ -  2.934(^-1)^   ... 


2i 


and 


,6.5)     ^=1.1^ 

TT 


3.812     -     1.763(^-1) 

+  1.235(4-1)^  -    .779(^-1)^ 


^  B-^(8.416)    J  Vgl-l 


V^  ^n  (/ITf  +/?:) 


(.00597)6^9 
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If  we  suppose  that  [o.h)   and  {'.'^)    remain  good  approximations 
when  the  parameter  ^.  which  we  assumed  to  be  small,  is  replaced  by 
e  =  1,  correspond'       rltical  flow,  then 


yp(i) 

"^    =  1  + 


^ 


f^^^^^]  -   ^  ^n  (/?+  1)  -  (.00597)0 

\  /2  +  1  /     TT^ 


=  .702  -  (.006)0  ,  0  <  0  <  1  . 

This  is  close  to  the  average  value  .71  noted  in  Section  4.   The 
acceptance  of  b  =  a[ .702  -  (.006)0]  as  an  approximation  to  the  brink 
depth  is  tantamount  to  the  assumption  that  X2(l)  is  close  to  zero. 
This  kind  of  assumption  is  implicit  in  some  of  the  analyses  which 
have  been  devised  for  the  calculation  of  b,  see  for  example,  Roy 
[4].   Consequently,  such  analyses  really  give  an  approximate  ordi- 
nate of  a  point  (a,b,)  on  the  upper  nappe  without  an  estimate  of 
the  presumably  small  abscissa  a.   According  to  our  second  order 
approximation,  b  =  a[ .702  -( .006)0]  is  the  depth  at 

X2(l)  =  f  [-  ^n  2  +  ^] 

IT 

=  -    .I4a  . 

In  general,  if  e  is  prescribed  we  can  invert  (6.4)  to  find 
i     such  that  x„(f,  )  =  0.   The  approximation  to  the  contraction 
coefficient  can  then  be  found  by  substituting  i^   in  (6.5)  in  order 
to  obtain  Y2^^    )/^-      An  estimate  of  the  brink  depth  can  also  be 
found  by  substltutinc  various  values  of  %   near  ^  =  1  in  (6.4)  and 
(6.5).   For  example,  take  e  =  .9  which  means  that  the  flow  is 


-!A 


To 


supercritical  to  the  extent  that  u  =  ^^-^ /ga  =  1.05  /ga.   For  this 
value  of  e  we  have 


^  -   (M).ll^-  (^---1-)' 


(6.6) 


and 


(6.7) 


X2(^)    1 


3   ••• 


-  .082[9.625  -  5.36^(e-l)  +  4.188(^-1)'^  ...]  J 


yo(?) 


a 


1  -  .009[3.8l2  -  1.763(^-1)  +  1.235(^-1)   .-.] 


+  .070  ^n(  ^^"^^^-M  -  ^-^^  ^n(/IT?:  +  /?:)-  (.005)0 


Computation  shows  that  for 


P  =  0  ,    e  =  .9 


we  have 


Xgd) 


061  , 


Yod) 


=  .725  -  (.005)9 


X2(.9) 


X2(.8) 


044  , 
.026  , 


yo(-9) 


yo(.8) 


=  .709  -  (.005)1 


=  .692  -  ( .005)0  ,         0  <  9  <  1 


A  comparison  of  these  figures,  and  others  we  have  obtained  in  this 
section,  with  the  values  noted  in  Section  4  Indicates  that  the 
approximations  (6.4)  and  (6.5)  are  rather  good  even  for  values  of  e 
close  to  1. 
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7.  Sec.i.j  Order  Approximation.   Equations  for  Shapes  of  Nappes 

We  have  already  seen  In  Section  6  how  the  approximate  shape 
of  the  upper  nappe  above  the  brink  can  be  found.   In  a  similar  way 
the  .     1  order  shape  of  the  lower  nappe  In  the  neighborhood  of 
the  brink  can  be  obtained  by  using  the  parametric  equations 


-^  =  Re  (CD^(U  +ea)^(^)  +  e%(a} 


(7.1) 


yo(^) 


=  1  +  Im  {a)Q(C)  +  ea)^(^)  +e^ai^{i)] 


where  i   is  restricted  to  -1  <  ^  <  0.   Collecting  previous  results 
for  o)  ,  a>,  and  oip  we  find  that  for  the  neighborhood  of  the  brink 
where  ^  =  -1  the  coordinates  of  the  lower  nappe  are  given  by 


(7.2) 


X2(0 


in  1^1  ^  £(p-fl)^n  \i\ 

IT  IT 


3(p+l)^^n  Ul 

27r 


+  e 


Mp+1) 


V.       TT- 


and 
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(7.3) 


yo(?) 


2n+l 


riP±lI  ^n    Ul    [^n    Ul   -4    in  2] 
27r'^ 


+    e 


^    (p+l)(8.4l6)    ^^      1-/1+^ 

TT  \  1  +vT+i; 


-  iit£^  -^  ta  (/IT?  +  /jIT) 

+  iPJ^  (2  in^  2  +  ^)  -  (.00597)(P+I)e 


where  0  is  a  real  number  such  that  0   <   9   <   1. 

The  shapes  of  the  nappes  downstream  corresponding  to  the 
neighborhood  of  ^  =  0  come  from  combining  (3-7),  (3.16)  and  (5.17) 
and  they  are  given  by  the  complex  form 


(7.4)      zgd) 


=   ai    -  ^  -  ^AJUA  _   aeiiP+il   [,^  ^  ^  ^  SttI    In  2] 


IT 


IT 


27r' 


(p+1)^  3  in  2        (p+l)^(8.4l6)    ,    (?+l)(8.4l6)  ^ 

2.5  .5 


+  ae 


2 ; 


6Tr^  271  5 


+ 


^^ii^    (27r^T   +   T^    -    14.424) 


6-5 


5  IT' 


V  ^  l(p+l)(8.4l6)T  _  l(p+l) 


TT 


pi  +  o(U|) 


y 


where   we  have   used   t   to  denote 
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T   =    In   ^    -    2    In   2    . 

The  upper  nappe  corresponds  to  positive  values  of  ^  while  the  lower 
nappe  corresponds  to  negative  values  of  ^.  Note  that  for  the  upper 
nappe  t  Is  real  and  negative,  while  for  the  lower  nappe 

T  =  iTT  +  j^-^-'- 

wlth  real  part  negative.   We  can  expect  (7-^)  to  be  a  good  approxi- 
mation if  e  is  small.   However,  the  approximation  is  not  good  when 
p  =  0  and  e  =  1.   Analysis  shows  that  the  breakdown  in  accuracy  for 
this  case  is  due  to  the  disappointing  circumstance  that  the  quad- 
ratio  terms  involving  t  cancel  out  of  (7-^)- 

It  is  interesting  to  see  what  {7-4)  implies  when  we  allow  a 
to  approach  zero  so  that  the  overfall  degenerates  into  a  sheet.   In 
order  to  maintain  a  constant  mass  flow  as  a  -►  0  the  velocity  u  must 
be  taken  so  that  au5  =  6h  is  constant.   Then,  if  we  set 

Pg  -  Pi  ax 

P  =   6h   '    ^  =  -  IRI 

we  find  that  as  a  -♦  0  the  real  and  imaginary  parts  of  (7-'+)  become 

X2(U  =  ua  -  p{pu  +  g)  -^ 

(7.5) 

y^ii)   =  -(pu  +  g)  ^ 


after  retaining  dominant  terms.   Now  it  can  be  verified  Independently 
of  what  appears  in  this  paper  that  in  the  first  stages  of  fall  the 
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shape  of  a  thin  overfall,  subject  to  gravitational  attraction  and 
different  constant  pressures  on  its  faces,  is  given  by 


X  =  M^  sin  Pt  -  $ 
P  P 

(7.6) 


=  M^  (cos  pt  -  1) 

p 


y ^ 


where  t  denotes  time.   In  expanded  form  we  have 

2 

X  -   ut  -  (g+pu)p'  -7—   .  .  . 

(7.7) 


-   t^ 

y  =  -(g  +pu) 


2"  •  ' 


This  agrees  with  (7.5). 
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